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Abstract 

It is shown that the use of a density dependent effective Pauh potential together with a nucleon- 
nucleon interaction potential plays a crucial role to reproduce not only the binding energies but 
also the matter root mean square radii of medium mass range spin-isospin saturated nuclei. This 
study is performed with a semiclassical Monte Carlo many-body simulation within the context 
of a simplified nucleon-nucleon interaction to focus on the effect of the genuine correlations due 
to the fermionic nature of nucleons. The procedure obtained is rather robust and it does not 
depend on the detailed features of the nucleon-nucleon interaction. For nuclei below saturation the 
density dependence may be represented in terms either of the nucleon number, A, or the associated 
Fermi momenta. When testing the simulation procedure for idealized "infinite" symmetric nuclear 
matter within the corresponding range of densities, it turns out that finite size effects affect the 
Pauli potential strength parametrization in systems up to about 120 particles while remaining 
approximately stable for larger systems. 
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I. INTRODUCTION 

We have recently demonstrated ^ that one can always expect to be able to reproduce 
the empirical binding energies for a set of nuclei by introducing a proper density dependent 
Pauli potential in terms of a single variable, the nucleon number A. Such a result is the 
consequence of a delicate counterbalance between kinetic and potential energetic contribu- 
tions. A deeper understanding of the nuclear saturation mechanism and reproduction of 
other observables for nuclei has been made possible due to the richness of available nuclear 
data. Since the pioneering shell model calculations [2^ and the so-called "Coester band" j^, 
many microscopic investigations have been devoted to study the properties of many-nucleon 
systems with a variety of approaches as i.e. relativistic mean fields J], quark-based mod- 
els js, ^ and chiral dynamics On the other hand, studies based on computer simulations 
with Monte Carlo techniques and/or (quantum) molecular dynamics, have proved them- 
selves to be v|er^ powerful to investigate many-nucleon systems such as subthreshold kaon 
production 



pasta phase 
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Olland multifragmentation [a, in heavy ion collisions, and the nuclear 
113|, llJ] . Most of these works, based on semiclassical simulations, treat 



nuclei as composite objects interacting through effective potentials containing the essential 
features of the system. Especially, the effects of the Pauli principle, the basic quantum 
feature for many-fermion systems, is simulated by means of an effective potential depending 
on the position and momentum of the interacting nucleons. This was firstly suggested by 



Wilets et al. 



15| in the context of classical many-body nuclear models. Later, approaches 



such as fermionic molecular dynamics (FMD) 16| and antisymmetrized molecular dynamics 
(AMD) 17|, include the Pauli principle in a quantum mechanical manner. However, it is 
still not practical to apply these methods to study heavy nuclei and nuclear matter, due 
to the need of heavy computational calculation, and unestablished (anti)periodic boundary 
conditions to simulate an infinitely large nuclear system 18|. 

In this work we go a step further from our previous study []]] by considering a refined 
effective spin-dependent NN interaction and analyzing other nuclear properties for medium 
mass nuclei and heavy nuclear systems. We will firstly deduce a corresponding density 
dependent effective Pauli potential reproducing the empirical nuclear binding energies. For 
this purpose, without loss of generality, we will focus on medium mass range nuclei with 
8 < A < 44. It turns out that the spin dependence has only a very moderate impact. 
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with the relevant contributions coming from the in-medium effects of the effective Pauh 
potential. In addition, we will also evaluate the root mean square radius (rms) of spin- 
isospin saturated nuclei and their density distributions. When the simulation is extended 
for idealized infinite symmetric nuclear matter, it turns out that finite size affects the Pauli 
potential strength up to about 120 particles while remaining almost stable for larger systems. 
The picture arising from the present semiclassical simulation of nuclei, though simplified, 
accounts for basic properties of these systems. In the medium mass range nuclei, the density 
dependence could be either described in terms of the nucleon number or, as we will show, a 
function of the associated nucleus Fermi momentum. This treatment relates to that of recent 

n 

approaches applied in chiral dynamics [7|. The present description can be considered, then, 
as an improvement of previous simulation work where neither antisymmetrization effects nor 
density dependence were considered. The derived effective Pauli potential does not depend 
on local densities as occurs, for example, in work using Skyrme models js]. Nonetheless, 
even if not local, the density dependence in the effective Pauli potential becomes crucial to 
achieve saturation. 

The structure of the paper is as follows. In section [Tll we describe the different nuclear 
potentials used in the simulations. First, we introduce 5*- wave NN interactions without spin 
nor isospin dependence. This makes easier to analyze the role of the Pauli potential and its 
nuclear density dependence. Next, we consider improved spin dependent NN interactions. 
Results will be presented in section IIIIl and summary and conclusions will be given in 
section [IVl 



II. NUCLEAR INTERACTION MODELS 



In this section we present the different simplified nuclear models which will be used in 
our simulations. We restrict ourselves to spin-isospin saturated, Z = N (even Z and A^) 
nuclei with mass number A in the range 8 < A < 44. In addition, nucleons are treated as 
classical, structureless particles. 
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A. Square-well potential 



First, we consider a NN interaction based on a square- well potential without any depen- 
dence on spin nor isospin. While the model is simple, it retains the basic features of the NN 
interaction, short-range repulsion and intermediate range attraction. The NN interaction 
potential between the i-th and j-th nucleons is given by, 



VNN{rij) = < 



Vcore, for < Tij < O, 

-Vo, for a < rij < a + 6, (1) 
0, for a -|- 6 < Vij, 



where Vij = |rj — rj| is their relative distance. This potential consists of a repulsive core 
of strength Vcore and width a, and an attractive well of depth Vq and width b. We will 
consider values Vcore = 100 MeV, Vq = 3 MeV, a = 2 fm and 6 = 2 fm. The core strength 
gives only a very moderate contribution to the binding energies (as easily deduced from the 
different values we may use [1,]), only preventing nucleons from occupying the inner region 
in configuration space. However the spatial parameters, a and b, are crucial to reproduce 
satisfactorily matter root mean square radii as will be explained later in this section. The 
corresponding Hamiltonian is given by, 

^ = E^+ E VNNiu,), (2) 

i=l i=ij>i 

where pj is the 3-momentum of i-th nucleon with mass mTv- From now on we will refer to 
this model as SW. 

Next, we consider a refined NN interaction including the Coulomb potential: 

Vcouiir,,) = ^(1/2 + r,)(l/2 + r,), (3) 

where Tj (tj) is the isospin third-component of i-th (j-th) nucleon (+1/2 for protons, —1/2 
for neutrons), and e is the proton electric charge. Then, the Hamiltonian is given by, 

^ = Eir^+ E [VNN{r^,) + Vcoul{r^,)]. (4) 
i=l 2"^iV i^ijyi 

We will refer to this model as SWC 

Since nucleons are fermions, they obey the Pauli principle. In the present treatment this 
is mimicked by introducing an effective potential, which prevents nucleons from occupying 



the same phase space volume when they have the same quantum numbers. This was first 
suggested by Wi.ets etM. Later, a Gaussian fonn of the Pau.i potential has been intrc 

duced by Dorso et al. [19|. In this way one can reproduce well kinetic energies per nucleon 
in a Fermi gas, although it fails to describe other basic features such as two-body correlation 



functions at low temperature [20]. Using this effective Pauli potential depending not only 
on the position but also on the momentum of the interacting fermions in a Hamiltonian 
approach, an effective nucleon mass which grows largely with increasing nucleon density 
arises. This results in very slow nucleons, and may lead to quasi-crystallized nuclear Fermi 



gases if not properly considered |2ll ]. 



In this work we use the Pauli potential form proposed by Dorso et al. [1 

- ^ 1 '^r,r/a,a, , (5) 

where pij = |pi — Pj| is the relative 3- momentum of the z-th and j-th nucleons, with Sr^rj 
(5o-jCTj) being the Kronecker's delta for the isospin (spin) third- component. We will differ- 
entiate two cases: (i) Vp, go and po are constant fixed to reproduce the empirical binding 
energy per nucleon (—E/A) of ^^O, 7.98 MeV J22I, and (ii) allow them to be density depen- 
dent. The models corresponding to (i) and (ii), will be referred to as SWCPo and SWCP, 
respectively. In addition, in both cases, spatial parameters in Eq. ([1]) take the values a = 2 
fm and 6 = 2 fm as they have been adjusted to reproduce the ^^O matter rms radius to 2.72 



fm 23|. The Hamiltonian for these models, without specifying the density dependence for 
the Pauli potential, is given by: 

H = '^^z-^+ [VNN{rij) + V{rij)coui + Vpauii{rij,pij)] . (6) 

i=l i=l,j>i 

To extract the density (Fermi momentum) dependence of Qq and Pq in the Pauli potential 



in Eq. ([5]), we use the Fermi momenta deduced in Ref. [2J| from the quasielastic electron 
scattering on several nuclei, and interpolate for the nuclei studied later. Note that with this 
parametrization of the Pauli potential the more general local density dependence has been 
parametrized in terms of one single parameter based on a simplified Fermi gas picture. 
In a nucleus, a typical nucleon sphere radius r may be given by. 



(7) 
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where p = 2p^p/37i'^ is the nucleon density and pp the Fermi momentum of nucleons in the 
nucleus. Then, the averaged inter-nucleon distance 2r may be estimated as (2r/v^5'o) — 1, 
where qo is an "effective range" of the Pauh potentiaL This leads to, 

;i(97r)i/3 

go (97r)V^ 

with the uncertainty principle for go and Pq to satisfy go Po — h. 

As mentioned above, we will first consider a model where the Pauli potential parame- 
ters, Vp, go, and po are fixed to reproduce the binding energy of ^^O, what gives rise to 
(Vp, go,po) = (41 MeV, 1.88 fm, 104.96 MeV/c). Secondly, we consider another parametriza- 
tion of the Pauli potential allowing a density (pp) dependence trying to reproduce the 
empirical binding energies for spin-isospin saturated nuclei, together with the relations in 
Eqs. dS]) and ©. This leads to. 



Vpauii{rij,pij,pp) = Vp{pp) exp 



/>i2 ^2 
ij Fij 



5r^rA.^^ (MeV). (10) 



H{Pf) 2pl{pp) 

It is worth noting that this approach has some similarities to what is used in other contexts, 
such as chiral dynamics [7]. This may imply that the variable pi;' in Eq. ( ITOl) . is appropriate 
to characterize the many-nucleon system. 



B. Spin dependent square- well potential 

The NN interaction potentials we have described so far may be more generally described 
by introducing spin dependence. Note that this also implies isospin dependence in order to 
satisfy the Pauli principle. We consider a simple potential for the spin parallel (TT, ii) and 
antiparallel (ti, it) pairs of nucleons: 

Vcore, for < Tij < a, 

(r.,) = -Vj^'''6.^^, - V,^'^'\l - for a < < a + b, (11) 

0, for a + 6 < Vij. 

We distinguish three different possibilities for the strength of each spin dependent term: 
(i) Model SDl, where the spin parallel potential is more attractive than the antiparallel, 
IK/^'^^I > iVo^'^^l; (ii) Model SD2, the opposite case, \V^^'^^\ < \V^^'^'^\; (iii) Model SD3, for 



TT I [ I T I I T I 

the case, | = |. We will consider values a = 1 fm and 6 = 2 fm, different from 

those used for the models described in Sect. Ill A[ In this way we introduce some dependence 
on the spatial parameter value a in Eq. ([T]). We have verified that parameter a (and also h) 
is mainly responsible for reproducing rms radii but has a very moderate impact on binding 
energies. We will mainly focus on the effect of the spin dependence on binding energies. The 
different balance between the parallel and antiparallel strengths will allow us to overcome 
the possible consequences of the simplification including only S'-wave interactions. Inclusion 
of also the Z)-wave interactions, which only couple to the ^Si channel (in the present case, 
e.g., the interaction between {p ]) — {n |) pair) would modify the balance between the 
attraction in the spin parallel and antiparallel cases. 

Finally, let us mention that although we have not considered a possible radial dependence 
on the interactions, in the case of binding energies will not have observable effects. The 
strength of Vq^^'^^ and V^^^'^^ for the SDl, SD2 and SD3 models are chosen rather arbitrarily 
but they all reproduce the empirical value, —E/A = 7.98 MeV of ^^O, together with the 
same Pauli potential in the SWCPo model. The Hamiltonian for these models is given by: 



H = J:7^+ E \VNTir,,) + Vp,^uir,„p,,) + Vcouiir.j)] . (12) 

i=l "^^N i = 



In Table H] we summarize the strength of Vg^^'^^ and Vg^^'^^ for all models. 

TABLE I: Strength (in MeV) of V^^'^^ and vj^'^'^ for the SDl, SD2 and SD3 models. 



Model 


Vcore 






SDl 


100 


5.8 


0.5 


SD2 


100 


0.5 


5 


SD3 


100 


3 


3 



III. SIMULATION RESULTS 

In this section we present the results obtained using Monte Carlo simulations of the 
nuclear systems where a typical value for low temperature, T = 1 MeV, is adopted. Initially, 
the nucleons are uniformly distributed inside a sphere of radius Rq, within a cubic box of 
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volume V = and impose L » Vij. Up to -Rq — 5 fm the binding energy per nucleon has 
rather stable values, although it tends to vanish rapidly for larger Rq. This reflects the fact 
that for large values of Rq some nucleons would not feel any attraction (see Eq. ([T])). Then, 



using the Metropolis algorithm 



251], the initial random seed is iterated until the energetic 



conflguration achieves the lowest possible value. In addition, we must check that a clustered 
conflguration of nucleons is indeed formed, a nucleus. Once the system is thermalized after 
performing a sufficiently large number (several thousands) of Monte Carlo sweeps, we take 
data to calculate the statistical average for thermodynamic quantities such as kinetic and 
potential energies and root mean square radii. 

First, we discuss the results of the SWCPo model introduced in Sect. Ill At based on the 
density independent Pauli potential fixed to reproduce the binding energy per nucleon and 
matter rms radius of ^^O. 

In Fig. [1] we show the SWCPo results for —E/A by the solid line compared to empirical 
values [22] (crosses). 

As can be seen, using the density independent Pauli potential, —E/A grows linearly as 
the nucleon number increases. Thus, the empirically observed saturation for —E/A cannot 
be achieved. 

Let us now discuss the results for the SWCP model, with the density dependent Pauli 
potential. Before presenting the results for —E/A., we show in Fig. [2] the pp for obtained 
to reproduce the empirical —E/A values for the spin-isospin saturated nuclei. The strength 
Vp{pf) increases as Fermi momentum (density) increases. The strength obtained for each 
nuclear species can be mapped as a function of the Fermi momentum extracted from the 
quasielastic electron scattering off nuclei 2J]. This behavior is analogous to the vector mean 
field in Hartree approximation in relativistic mean field models j^]. 

We show in Fig. [3] the results for —E/A obtained in the SWCP model with the use of the 
density dependent Pauli potential, Eq. ( fTOl) . The statistical uncertainty of these calculated 
values is at most 5 %. 

In Fig. Hlwe show the matter rms radius (boxes) as a function of mass number. A, for a 
set of nuclei as calculated with the density dependent fit in Eq. ffTOj) . We can see that they 
approximately follow the liquid drop model results (dashed line) given by the R = tqA^^^ 
with tq ^ 1.1 fm. 

In Figs. O and Owe show the density distribution calculated for an ^^O nucleus and a 
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FIG. 1: Binding energy per nucleon, —E/A, calculated in the SWCPo model. Empirical data, 
shown by crosses, are taken from Ref. 



22( 1 for even nuclei. We have calculated nuclei with A = 
8, 12, 16, ...,44 and the solid line goes through all calculated points, being only valid for these even 
nuclei. This will also apply to the other figures along the text. 



^°Ne nucleus, respectively. We can see that most of the mass is located in the nuclear central 
region. The decrease of particle population at small distances is due to the strong repulsive 
core of the interaction used. We have fit this distribution to a typical Woods-Saxon form 
p(r) = ^^^(/L^)/^ obtaining a set of values for ^^O, (3 = 0.066 ± 0.002/m-3, a = 4.0 ± 0.5/m 
and 7 = 0.41 ± 0.45/m and for ^^Ne we get (3 = 0.068 ± 0.007/m-^ a = 3.97 ± 0.20/m and 
7 = 0.1 ± 0.5/m which are in reasonable agreement with those of mean field calculations {4]. 
Note that the value a indicates where the central density remains approximately constant, 
while the matter rms radius gets contribution also from the more extended spatial region 
of nucleons in the nucleus. It is worth noting that the values of the parameter a and b are 
mainly responsible for reproducing the rms radii in this approach. 

Next, we present in Fig. [7] the comparison for the binding energy per nucleon obtained 
with the different models described in Sect. Ill A[ The sohd (long-dashed) [short-dashed] 
hue shows the SW (SWC) [SWCP] model. Both, the SW and SWC models, cannot achieve 
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FIG. 2: Pauli potential strength Vp as a function of the Fermi momentum, pp- 
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FIG. 3: Binding energy per nucleon for spin-isospin saturated nuclei for the SWCP model. 
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FIG. 4: Matter root mean square radius calculated as a function of mass number A for several 
nuclei using the SWCP model. See details in the text. Liquid drop model calculation is shown for 
comparison. 



the saturation for —E/A. The Coulomb potential can negligibly help to achieve —E/A 
saturation as expected, although it is important for some properties of finite nuclei ^. 

In order to understand how the energy saturation is obtained in the SWCP model, we 
draw separately in Fig. [8] the contribution of the kinetic energy per nucleon K/A (dashed 
line), and potential energy per nucleon V/A (short-dashed line) together with E/A (solid 
line). We see the kinetic energy per nucleon largely deviates from the classical value, 
K/A = (3/2)fcsT, (i.e. K/A = 1.5 MeV in the present case with T = 1 MeV) due to 
the contribution of the density dependent Pauli potential. Thus, this shows that the kinetic 
and canonical momenta are completely different as the nucleon density increases [21.1]. A 
balance between the large potential and kinetic energy contribution leads to —E/A satura- 
tion. This cancellation arises as a consequence of the increasing kinetic contribution of Pauli 
correlations [l|. This may be analogous to the mechanism of saturation in relativistic mean 
field models which involves the cancellation between the attractive scalar and the repulsive 
vector potentials {4]. 
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FIG. 5: ^^O density distribution histogram as given by the SWCP model. Dashed Hne shows the 
Woods-Saxon fit for this nucleus. See details in the text. 
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FIG. 6: Same as Fig. [5] for ^"Ne. 
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Let us now study the effect of spin dependent potentials, using the SDl, SD2 and SD3 
models described in Sect. IIIBI In this case we will work with the density independent Pauli 
potential described in Sect. Ill Al in order to see the effect of spin-dependent forces in the 
energy saturation. The results are shown in Fig. [9l The solid (long-dashed) [short-dashed] 
line stands for the results of the SDl (SD2) [SD3] model, together with the empirical data 
(crosses) 22|. We can see again the —E/A saturation cannot be obtained even using an 
improved spin dependent NN potential. 

In Fig. [To] we show energy per nucleon for the ideal " infinite" symmetric nuclear matter 
(SNM) case at T = MeV (solid line). An incompressibility value of ~ 230 MeV has been 
assumed and po = 0.165/m~'^. Correspondingly the range in densities or Fermi momentum 
can be extracted as done in ref . [24|] as p = for the set of symmetric spin saturated nuclei 
8 < A < 44 (dashed line). The lower and upper bounds in the density interval correspond 
to p/po = 0.51 and p/po = 0.8 respectively. 

We now compare some results from heavier nuclear systems than the medium mass range 
nuclei considered in the previous sections in this work. In Fig.[TT]we show the Pauli potential 
strength, Vp as a function of the particle number. A, for an idealized case of SNM, through 
the simulation for densities p/po = 0.58 (solid line) and p/po = 0.7 (dashed line). We can 
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FIG. 8: Kinetic {K/A) and potential (V/A) energy per nucleon contributions to E/A. 
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FIG. 10: Energy per particle in SNM at T = and range in densities corresponding to the Fermi 
momentum of nuclei considered in this work. The lower bound in the density interval corresponds 
to p/po = 0-51 and the upper to p/ po = 0.8. 



see that finite size effects (particle number) are important up to values about A ^ 120 where 
the Pauli potential strength stabilize to an approximately constant value. The larger values 
for small A do not match those of finite nuclei since in this case nucleons are artificially 
constrained in a cube of length L = {A/ pY^^. We use periodic boundary conditions as usual 
to minimize boundary effects near the walls of the simulation box. It is worth mentioning 
that as we consider heavier nuclei pj] the density dependence parametrized by the Fermi 
momentum must be replaced by the system size, A, since one can map out Pauli potential 
strengths in a clearly defined way. In Fig. [12] we show a thermalized configuration of a 
simulation of a SNM plasma at T = 1 MeV at a density p/po = 0.7. In this homogeneous 
system of A = 200 particles, protons are shown as small dots while neutrons correspond to 
larger dots. 
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FIG. 11: Pauli potential strength as a function of the simulation system particle number for SNM 
for densities p/po = 0.58 and p/po = 0.7. 




FIG. 12: Simulation box for a system of SNM consisting of 200 particles at a density p/po = 0.7 
and r = 1 MeV. Protons and neutrons are depicted with small and larger dots respectively. 
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IV. SUMMARY AND CONCLUSIONS 



Using many-body simulations with Monte Carlo techniques we have pursued the study of 
an effective Pauli potential in the nuclear binding energy and matter rms radius considering 
density dependence. The density dependent effective Pauli potential partially simulates the 
effect of genuine fermionic correlations on semiclassical descriptions. For medium mass range 
nuclei with 8 < A < 44 the nucleon number dependence into the effective Pauli potential 
can be replaced in terms of the Fermi momentum of the constituent nucleons. Even if spin 
dependent effects are explicitly included in the NN potential, an appropriate parametrization 
of the density dependence can still be derived in terms of the Fermi momentum. The density 
dependence in the effective Pauli potential is crucial to reproduce nuclear binding energy 
saturation, wich is achieved by a balance of density dependent attractive and repulsive 
contributions in an analogous way to the relativistic mean field approach. Although for 
medium mass nuclei the density dependence can be described in terms of the Fermi momenta, 
the treatment of heavier nuclei requires to consider the density dependence in terms of the 
nucleon number, due to saturation of the Fermi momenta around 270 MeV/c. 

The use of spin-isospin dependent NN potentials turns out to give only a moderate 
improvement for the binding energy per nucleon, but it is unable to achieve energy saturation 
unless a density dependent Pauli potential is used. Our results show that, provided a different 
set of values for the spatial parameters, a and b or potential well strength, Vq, in the NN 
potential, app dependent Pauli potential strength can always be parametrized to reproduce 
the empirical nuclear binding energies. The procedure presented in this work is robust. 
Then, each set of nuclei modelled by a given NN potential and calibrated by the empirical 
binding energies, can be tested by studying further properties of nuclei like, for example, 
when introducing some hypcronic content. 

We have shown that, in the density range close to the saturation of nuclear binding 
energies, and for an idealized SNM system the Pauli potential strength depends on the size 
of the system and stabilizes to an approximately constant value around A 120. The 
matter rms radii of medium mass range nuclei can be satisfactorily reproduced by tuning 
the spatial parameters of the NN interaction to a reference nucleus (in this work taken to be 
^^O). A close agreement with the typical sizes of nuclei expected from a liquid drop model 
is obtained. 
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To summarize, the present approach may provide a practical efficient way to model a set 
of spin-isospin symmetric nuclei calibrated by the empirical binding energies and matter rms 
radii. It also can be regarded as a preliminary step for studying further properties of sym- 
metric many-nucleon systems in semiclassical simulations. It is necessary to provide further 
insight into semiclassical effective potential models to study wether asymmetric systems can 
be treated in a similar fashion. 
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